Noise-driven quantum criticality 



O 
(N 

o 

Q 



Oh 



> 
o 



X 



J. Eisert^ and T. Prosen^'^ 
^ Institute of Physics and Astronomy, University of Potsdam, 14476 Potsdam, Germany and 
^Department of Physics, FMF, University of Ljubljana, 1000 Ljubljana, Slovenia 

We discuss a notion of quantum critical exponents in open quantum many-body systems driven by quantum 
noise. We show that in translationally invariant quantum lattice models undergoing quasi-local Markovian 
dissipative processes, mixed states emerge as stationary points that show scaling laws for the divergence of 
correlation lengths giving rise to well-defined critical exponents. The main new technical tool developed here is 
a complete description of steady states of free bosonic or fermionic translationally invariant systems driven by 
quantum noise: This approach allows to express all correlation properties in terms of a symbol, paralleling the 
Fisher-Hartwig theory used for ground state properties of free models. We discuss critical exponents arising in 
bosonic and fermionic models. Finally, we relate the findings to recent work on dissipative preparation of pure 
dark and matrix product states by Markovian noisy processes and sketch further perspectives. 



Quantum phase transitions (QPT) and quantum criticality 
are among the most remarkable emergent features and prop- 
erties of large quantum quantum lattice systems. Even if indi- 
vidual constituents directly interact only to their immediate 
neighbors, characteristic length scales associated with two- 
body correlation functions diverge at critical points. Critical 
phenomena are responsible for a sizable number of effects in 
quantum many -body physics llllll]- Specifically, if the Hamil- 
tonian is of the form 



(1) 



where ft, is a Hamiltonian term acting on a small number of 
sites only — usually nearest neighbors — and Tj is the shift to 
the lattice site j, and V is another Hamiltonian of this form, 
then a system is critical at g = gc if the ground state energy of 
H = Hq + gV is non-analytical at gc- In second order QPT 
this is accompanied with a divergence of coiTelation lengths 
giving rise to the familiar notion of a critical point in quantum 
many-body systems. 

In this work, we quantitatively explore the possibility of 
scaling laws of correlation lengths of quantum many-body 
systems driven by quantum noise. The Hamiltonian will re- 
main entirely unchanged, merely we add Markovian noise that 
acts quasi-locally, so either on single sites alone or on a fi- 
nite number of neighbors. We will see that despite this lo- 
cal action of the dissipative dynamics, mixed states arise as 
steady states that show characteristic divergences of correla- 
tion lengths, even allowing for the meaningful definition of 
critical exponents. The observation that local noise alone can 
be used to prepare pure quantum many-body states exhibit- 
inglong-range order was made in the seminal work of Refs. 
ISl-lg]. Here, we go beyond that insight by developing a theory 
of critical exponents of (generic ally non-pure) non-equilibrum 
steady states of local dissipative systems. The main new tech- 
nical tool will be a description of free translationally invari- 
ant bosonic and fermionic open quantum systems driven by 
Markovian quantum noise — paralleling the central role free 
models play in the common theory of QPT. This approach al- 
lows to express correlation functions of canonical coordinates 
in terms of a symbol reflecting a momentum space descrip- 



tion, reminding of the so powerful Fisher-Hartwig formalism 
ItIi for ground state properties. 

Noise-driven critical exponents. As far as the Hamiltonians 
are concerned, we study the — in the theory of criticality — 
usual local quantum many-body systems on a lattice. The 
Hamiltonian can hence be written as in Eq. ([T]i; this is nothing 
but the familiar general translationally invariant local lattice 
model. In addition to the local Hamiltonian, we allow for 
noisy processes merely locally affecting the quantum many 
body system. This noisy process is supposed to be Marko- 
vian. Hence, each site j is associated with a generator of form 

^aiP) = lY.i^^^(9)pL,ig)^ - {L^{g)^L^ig),p}) , 

where g is a parameter in the many-body model. The Lind- 
blad operators {Lfj,{g)} in turn act merely on finitely maiiy 
sites, typically on-site or nearest neighbors. It is key to this 
work that noisy processes with long-range interactions are 
excluded — and still we will find divergent correlation func- 
tions. The equations of motion of the state are given by 



-p=-i^[r,(ft),p]+^r,(/:,(p)) 



(2) 



where the summations (j) run over the d-dimensional lattice 
Z'^. For local observables {Oj}, where Oj is supported on site 
j only — again as usual — one considers correlation functions 

f{j,k)^\{0,Ok)-{0,){Ok)\ (3) 
and defines the correlation length as the length scale 

C = - lim logf{j,k)/\j-k\. (4) 

In sharp contrast to Hamiltonian standard descriptions, we 
now put into the center of attention the dependence of steady 
states in terms of the noise only, while keeping the Hamilto- 
nian unchanged. As the main focus of this work, we will see 
that it makes sense to define a critical value gc of this param- 
eter, such that when approaching 



9^ 9c 



(5) 
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the correlation length will diverge. In fact, one can define a 
noise-driven critical exponent A, by means of the asymptotic 
behavior of this correlation length 



A|.g-.gc|' 



(6) 



for some A > 0. Note that while this parameter is defined as 
the critical exponent for critical phenomena in the usual sense, 
all divergences of correlation lengths are entirely driven by 
quantum noise in this open quantum system. 

Free driven models. Centre stage will be taken by free mod- 
els, so models where the Hamiltonian h is quadratic in bosonic 
or fermionic operators and L^{g) are linear Bosonic and 
fermionic systems will be treated in an as parallel as possible 
fashion fl. Bosonic operators {hj} associated with site j are 
accompanied by commutation relations, whereas fermionic 
operators {fj} anti-commute. 



(7) 



Suppose for a moment we ignore the translational invariance 
and consider a general quasi-free many body system defined 
on n sites, where the limit n oo will be considered later 
The canonical variables can be grouped into a 2?! dimensional 
vector of positions and momenta for bosons u = (u^j ; ly = 
1,2; 1 < j < n), where uij = bj + 5^, = \{bj — 6p 
and analogously of Majorana operators w for fermions with 
wij = fj + /| and W2.j ~ — /|). This convention, 
different from the usual one by a factor of \/2, simplifies some 
of the later expressions. 

We (i) assume the Hamiltonian ([T]i to be defined by a 



p2nx2n 



quadratic form, H = u H^u, where ~ G 
for bosons, and H ~ w'^HfW where Hf 
purely imaginary, matrix for fermions. (ii) The Lindblad op- 
erators = /Ju (L^ = ^Jw;) are assumed to be linear in 
the coordinates and can be written in terms of 2?! dimensional 
vectors l^, which can be grouped together into a Hermitian 
bath matrix 



(8) 



M E {£2nx2n^ having a symmetric or anti-symmetric real 
or imaginary part, = {M + M)/2 = and Mi 
— i(A/ — M)/2 ~ — M^, respectively. The elementary two- 
point correlations are gathered in a Hermitian covariance ma- 
trix, being a real, symmetric matrix for bosons 7b G ]g2nx2n^ 

(7b)j,fc = ^ ti-p{ujUk + UkUj) = Rctr(/?UjUfc), 



and a real anti-symmetric matrix for fermions 7f G 



p2Tix2n 



(7f)j\fc = i^^'^PiwjWk - WkWj) = -ImtripWjWk). 

Obviously, if one adds terms linear in bosonic operators to the 
Hamiltonian (fermionic first moments are always vanishing), 
then one should subtract the possibly non-vanishing mean 



values in the definition of the covariances. Defining a sym- 
plectic unit matrix cr = iay (E) In, the bosonic covariances 
satisfy the uncertainty relations in the form 7b + icr > 0, 
whereas an analogous relation for the fermionic covariances 
reads ■yf + l2„ > 0. 

From the equations of motion of the canonical coordinates 
or the Majorana operators in the Heisenberg picture, one finds 
after a tedious but straightforward computation that the co- 
variance matrix satisfies a closed set of equations of motion 



dt 



(9) 



taking the same form for bosons and fermions. The specific 
matrices governing this equation of motion are given by 

Xb = a{2Hb + 2Mi), Fb = M^a = Y^^ (10) 



for bosons and 



Xi = -2iH{ + 2M„ Yi = 4i\/i = -y/ 



(11) 



for fermions, where Xh.,Xi,Yh,Yi G R^nxin^ ^oie, an im- 
portant distinction between bosonic and fermionic quasi-free 
semigroups: Namely the spectrum of fermionic Xf (which de- 
termines all the decay rates of Liouvillean relaxation Igi lToll ) 
can only lie on the positive real-side of the complex plane 
(since A/j- > 0, see Lemma 2.3 of Ref. ifioll '). whereas in the 
bosonic case the spectrum of X\j dTol ) (which again determines 
all the Liouvillean decay rates [1 lj|) is not constrained, so the 
Liouvillean evolution becomes unstable — describing indefi- 
nite pumping-in of energy — when some of the eigenvalues of 
X\, attain negative real part. It is also easy to identifypure 
steady states, corresponding to the dark states of Ref s. ||3l|5|], 
but here applied to the free setting: They have covariance ma- 
trices for which all eigenvalues of (cr7b)^ or of (7f)^ are all 



equal to 1 in the bosonic and fermionic case, respectively Ill2ll . 

Translationally invariant systems and symbols. Next, we 
consider explicit translational invariance on j G Z^, so n = 
L"^, and subsequently take the limit L — > oo. In case of finite 
systems, the yet to be defined symbol will be defined over 
a discrete quasi-momentum space and the analysis does not 
apply: Hence, in all that follows, divergences of correlation 
lengths manifest genuine quantum many-body effects. Due to 
translational invariance we may write 

iH,,)(u,j),{iy'.j') = {hrjU v,v' 1,2, 



and similarly for the Lindblad operators L^j introducing vec- 
tors on Z^. Then we may define the Hamiltonian 
and the Lindbladian symbols, as, respectively, 2x2 matrix 
valued and 2-dimensional vector valued functions on T*^ = 



9 ^ 



K{3)e-'^-\ 



Note that here denotes only those Lindblad vectors which 
directly couple to the site = of the lattice, all the other 
Lindblad vectors are obtained by the translations (as in Eq. 
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dU). By the convolution theorem, the bath matrix M of Eq. 
(O is then also circulant, with a symbol 



which will be, for convenience, decoupled into a real and 
imaginary part TOr(<y9) = {m{if) + {—ip))/ 2 and fhi{tp) = 



-i{m{ip) 



{—ip))/2, respectively. 



Furthermore, if the stationary point of the dynamical semi- 
group (O is unique, it is also translationally invariant: This 
follows immediately from a group twirl with shifts tj form- 
ing a unitary representation: The steady state covariance ma- 
trix 7,, is also (block) circulant. Then the condition for the 
fixed point of Q, namely the Lyapunov-Sylvester equation 



1 1311 X^jjj + jjjX = Y^i maps to a simple 2x2 matrix equa- 



tion for the covariance symbol 



exp(-iv5-j), (12) 



namely 



(13) 



4crJmiCr2 



where for bosons ib — a2{2h\, + 2mi),yb 
(where a2 — icTy), and for fermions if = ~2ihf + 2rh^, yi — 
Ami in direct consequence of (IIOII 11 1. We note that Eq. (fT2T i is 
a simple 2x2 Sylvester matrix equation (in fact a linear sys- 
tem of 4 X 4 equations for matrix elements of 7^) for a given 
fixed quasi-momentum ip, and can be trivially solved. The 
real-space two-point correlation function (|3]l is then obtain by 
inverse Fourier transforming the covariance symbol 



1 



{2TrY 



d'^ip%{ip) exp(iv3 • r). 



(14) 



One-dimensional models. Clearly, for systems with short- 
range or finite-range interaction, the symbols h^i and to, and 



hence i^, are holomorphic on 



therefore the solution 



%{<f) of (fT4l) can have at most simple pole singularities, 
(p* G C^, determined by either of the following conditions 

f3,{p*)+M-^n^0. or /3i(^*)+/?2(-^*) = 

where f3i_2{'p) denote the two eigenvalues of Xn{p). For the 
rest of the discussion we focus on one dimension d = 1. Then, 
if ip* is an isolated pole singularity with minimal |Im p* |, then 
it takes straightforward complex analysis to show that 3C > 
0, such that 



lim 



00 exp(— |Im(/3*||r|) 



C, 



IW|.(15) 



This comprises the central result of our Letter The criticality 
of the system is then signaled by the closest pole p* approach- 
ing the real circle T^. Clearly, due to continuity of x{(p), this 
means that for a Bloch quasi-momentum ip' = Reip*, the 
eigenvalue of x{p)), whose real part l/r determines the relax- 
ation time of the corresponding Liouvillean normal mode, can 



be estimated as l/r < C"|Im<y9*| where constant C may de- 
pend on system parameters but not on (p* . Thus the criticality 
(|6]l also implies a critical slowing down r > C"^. Note that 
this way of finding critical points is different from the usual 
approach in Hamiltonian problems, constituting a difference 
in how dissipative and non-dissipative systems are treated. 
Dynamical instabilities may occur and may be detected in this 
framework, but the divergent correlation lenghts always re- 
flects a true quantum many -body effect and can not be viewed 
as local dynamical instabilities. 

For the fermionic case, it can also be shown that in case 
the pole arrives right down to the real torus T'', i.e. Inup* = 
0, then the solution of Sylvester equation there has to vanish 
7f(<^' = p)*) = (a consequence of Lemma 2.5 of Ref. llOll ). 
and hence algebraic decay of correlations is not achievable. 

Example 1: Driven criticality for fermions. We conclude 
the discussion of the quasi-free case by providing explicit ex- 
amples. An interesting example in the fermionic case in one 
dimension {d = 1) is provided by the fermionic model cor- 
responding to a XY spin-1/2 chain with exchange couplings 
Jx = ^(1 + r), Jy = i(l — r) and transverse magnetic field 
B, by virtue of the Jordan- Wigner transformation (see, e.g., 
Ref. lfl4ll for the exact description of the model). The Hamil- 
tonian symbol is hf{ip) — ^{B — cosip)(Ty + E-(sin </j)(Tx- 
We first discuss the case of a spatially incoherent, on-site 
noise, where a single Lindblad operator per site is considered 
Lj{g) = e{f, + /]) + d(/, - /j)e'f, where we find 



r'=Im^*~e2((|S|-l) 



4 • 2 

e sm 



.9). 



This means that we may have criticality as 5 — > = 
or gc = TT, only if the (noise-free) Hamiltonian XY model 
is already critical, namely if |-B| = 1. Therefore, on-site 
noise cannot induce criticality. However, things change drasti- 
cally when we consider two-site spatially coherent noise. For 
simplicity we again take a single Lindblad operator per site 
Ljig) = + /]) + £(/,+! + /j+i)e'». This term does 
not violate the parity of fermion number superselection rule 
lITsIl and defines valid and physically meaningful noise. In 
this case, our analysis yields a very simple expression for the 
covariance symbol 



7f((p) 



sm g sm ip 



2(1 + cos 5 cos (/s) 



which, quite astonishingly, does not depend on any of the 
hamiltonian parameters (F, B) at all! Now, we find a clear, 
noise induced criticality for gc = ot gc = tt, with the in- 
verse correlation length = arcosh(l/ cos g), which close 
to the critical point can be estimated as ^ ^ \g—gc\~^^^ giving 
the critical exponent A = 1/2, same as in the boundary-driven 
XY model ll4ll . and also as for ground states of free bosons. 

Example 2: Driven criticality for bosons. As a second ex- 
ample we consider, for t,v E M, a simple nearest neighbor 
hopping translationally invariant bosonic chain, 

H = + b]+,b, -{v + 2)b]b,), 

3 



4 



with Hamiltonian symbol h\^{(p) = t{cos(p — v), and on-site 
noise with a single Lindblad operator per site Lj{g) — e{bj + 
fej) + ei{bj — 6j)e'f in analogy to the previous fermionic case. 
The eigenvalues of X^, /3i.2 — sing ± 4it|w — cosiy9|^/^ 
signal instability of the Liouvillean fixed point (steady state) 
for g £ (— TT, 0) where Re /3i_2 < 0. Following previous para- 
graphs, the covariance symbol reads 



7b ('^) 



^ 1 .2te (^) — cos cos p 
sin g ' 



sin g cos g 



1 

sing 



e sm g cos g ^ 
2te {v— cos If) cos g 



z{ip) = 4i^(u — cosifY + £'*sin^g. The pole singularity 
is determined by a zero of z{lp*) ~ Q, Lp* ~ arccos(u ± 
\esmg/2t). Two regimes emerge: (i) If one has an optical 
gap \v\ > 1, Im ip* is bounded from below independently of 
g and no noise-induced criticality occurs, (ii) In the acoustic 
regime \v\ < 1, yet, we may expand Imip* in sin g and obtain 



(e/(2^(l-^;2)l/2))|sing|, 



giving the critical points gc = 0, tt and the critical expo- 
nent A = 1. We note, however, that noise induced criticality 
of a bosonic system corresponds exactly with the instability 
thresholds discussed above. A similar analysis can be imple- 
mented also for two-site coherent noise. 

Remarks on the relationship with matrix-product states as 
dark states. The works fsl [st] consider pure states as dark 
states of dissipative processes. Indeed, following the ideas 
of Ref. |5,], quasi-local dissipative Markovian processes can 
be constructed with matrix product states llol being the 
unique dark states. This insight can be combined with one 
of Ref. 1171], where parent Hamiltonians are constructed hav- 
ing unique ground states, and which undergo a quantum phase 
transition of any order when parameters of the Hamiltonian 
are being altered. E.g., Ref. lIlTIl considers H = Tj(h), 
where, with S being the vector of spin-1 matrices, 

/i = (2 + ® 5) + 2(5 ® Sf + 2(4 - g''){S^ ® if 

+ [g + 2f{S, ® S,f + g{g + 2){S, ®S,,S® S), 

with unique ground state for g gc = 0- Following the 
construction of Ref. one can then identify quasi-local dis- 
sipative processes - in just the same sense as above - giving 
rise to critical exponent A = 1 flsll. 

Steady states and entanglement area laws. We finally dis- 
cuss the entanglement structure of steady states driven by dis- 
sipation and identify a connection to entanglement area laws. 
This connection is particularly manifest in the case of free 
bosonic systems. Since their steady states are usually mixed, 
a genuine entanglement measure, and not the entropy of a 
reduced state, has to be applied. We bound these entangle- 
ment measures of a region A distinguished from the rest of 
a one-dimensional chain L, such as the distillable entangle- 
ment, from above by the logarithmic negativity lll9n . defined 



for a state p as Em{p) = log2 \\p^ where denotes the 
partial transpose of p with respect to the degrees of freedom 
of A. The state being Gaussian, the covariance matrix of p^ 



is given by 7^ = P'jhP, where P — 1„ ® Q with Q be- 
ing the diagonal matrix whose main diagonal elements are — 1 
when a degree of freedom is associated with A and 1 oth- 
erwise. The logarithmic negativity can in turn be computed 
from the symplectic eigenvalues {Xj} of 7^^^, given by the 
singly counted positive square roots of the eigenvalues of the 
matrix K = (crY^^^Eio'lEi—o'Y^^- Using a bound remind- 



ing of one established in Ref. 



, one finds 



log2 = ^log2max(l, A7I) < J2 IV' - 1| 



|A'-^/'-l2J|i<l|A-i/2 



where the norm || adds all absolute values of all matrix en- 
tries. From this, one can see that whenever entries of 7b (j^) 
are exponentially decaying away from the main diagonal — as 
they always do in the models considered here — one encoun- 
ters an upper bound that does not depend on the number of 
sites of A. Hence all steady states of the considered one- 
dimensional bosonic models satisfy an entanglement area law. 

Summary and outlook. In this work, we introduced the 
concept of critical exponents of quantum many-body systems 
driven by quantum noise. We identified scaling laws remind- 
ing of the theory of criticaUty for ground states of local quan- 
tum many-body models. At the heart of our discussion is a 
theory of driven free models — providing an important simple 
class of models for which all features of criticality can be es- 
tablished analytically, reminding of the theory of ground state 
criticality for free models. It is the hope that this work con- 
tributes to paving the ground for a systematic study of prepa- 
ration of complex states by dissipation and quantum noise. 

Acknowledgements. We warmly thank S. Diehl for inspir- 
ing discussions and feedback on the manuscript and the EU 
(Qessence, Compas, Minos), the EURYI scheme, and the 
Alexander von Humboldt Foundation for support. 



[1] S. Sachdev, Quantum phase transitions (Cambridge University 

Press, Cambridge, 1999). 
[2] M. Vojta, Rep. Prog. Phys. 66, 2069 (2003). 
[3] S. Diehl, A. Micheli, A. Kantian, B. Kraus, H. P Buchler, and 

P Zoller, Nat. Phys. 4, 878 (2008). 
[4] F. Verstraete, M. M. Wolf, and J. I. Cirac, Nat. Phys. 5, 633 

(2009). 

[5] B. Kraus, H. P Buchler, S. Diehl, A. Kantian, M. Micheli, and 

P Zoller, Phys. Rev. A 78, 042307 (2008). 
[6] S. Diehl, A. Tomadin, A. Micheli, R. Fazio, and P Zoller, Phys. 

Rev. Lett. 105, 015702 (2010); S. Diehl, W. Yi, A. J. Daley, and 

P Zoller, ibid. 105, 227001 (2010); A. Tomadin, S. Diehl, and 

P Zoller, arXiv: 101 1.3207. 
[7] A. R. Its and V. E. Korepin, J. Stat. Phys. 137, 1014 (2009); 

B.-Q. Jin and V. E. Korepin, ibid. 116, 79 (2004); J. Eisert and 

M. Cramer, Phys. Rev. A 72, 042112 (2005); P Calabrese and 

F H. L. Essler, arXiv: 1006.3420. 
[8] M. Cramer, J. Eisert, and M. B. Plenio, Phys. Rev. Lett. 98, 

220603 (2007). 
[9] T. Prosen, New J. Phys. 10, 043026 (2008). 



5 



[10] T. Prosen, J. Stat. Mech. P07020 (2010). 

[11] T. Prosen and T. H. Seligman, J. Phys. A 43, 392004 (2010). 

[12] To provide a complete classification of free dark states consti- 
tutes an interesting further perspective. 

[13] R. H. Bartels and G. W. Stewart, Comm. ACM 15, 820 (1972). 

[14] T. Prosen and 1. Pizorn, Phys. Rev. Lett. 101, 105701 (2008). 

[15] Only when viewed as a spin model, the noise becomes non- 
local as Jordan- Wigner strings arise. 

[16] M. Fannes, B. Nachtergaele, and R. F. Werner, Commun. Math. 
Phys. 144, 443 (1992); D. Perez-Garcia, F. Verstraete, M. M. 
Wolf, and J. 1. Cirac, Quant. Inf. Comp. 7, 401 (2007). 



[17] M. M. Wolf, G. Ortiz, F. Verstraete, and J. I. Cirac, Phys. Rev. 

Lett. 97, 110403 (2007). 
[18] This argument as such does not provide a bound to the spectral 

gap dictating the rate of convergence to this stationary point. 
[19] J. Eisert, PhD thesis (Potsdam, February 2001); G. Vidal and R. 

F Werner, Phys. Rev. A 65, 032314 (2002); M. B. Plenio, Phys. 

Rev. Lett. 95, 090503 (2005). 
[20] J. Eisert, M. Cramer, and M. B. Plenio, Rev. Mod. Phys. 82, 

277 (2010); M. B. Plenio, J. Eisert, J. Dreissig, and M. Cramer, 

Phys. Rev. Lett. 94, 060503 (2005). 



